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Abstract 



The Hamiltonian structure of the monodromy preserving deformation equations of Jimbo et al [JMMS] 
is explained in terms of parameter dependent pairs of moment maps from a symplectic vector space to the dual 



spaces of two different loop algebras. The nonautonomous Hamiltonian systems generating the deformations 
are obtained by pulling back spectral invariants on Poisson subspaces consisting of elements that are rational in 
the loop parameter and identifying the deformation parameters with those determining the moment maps. This 
construction is shown to lead to "dual" pairs of matrix differential operators whose monodromy is preserved 
1 under the same family of deformations. As illustrative examples, involving discrete and continuous reductions, a 

On ■ higher rank generalization of the Hamiltonian equations governing the correlation functions for an impenetrable 

Bose gas is obtained, as well as dual pairs of isomonodromy representations for the equations of the Painleve 
O !* transcendents Py and Pyi- 

1. Monodromy Preserving Hamiltonian Systems 

•rH , 

The following integrable Pfaffian system was studied by Jimbo, Miwa, Mori and Sato in 
[JMMS]: 

n 

dN, = - J2Wi, Njjdlogiai - a,) - [N t , d{a t Y) + Q}. (1.1) 

Here {Ni(ai, . . . , a n , y±, . . . , y r )}i=i,... n is a set of r x r matrix functions of n + r (real or 
complex) variables {ai,y a } <=i,...n , Y is the diagonal r x r matrix Y = diag(yi, . . . ,y r ) and 

a = 1 , . . . r 

the matrix differential form G is defined by: 

n 

e a6 = (1 - S ab )Q2Ni) ab dlog(y a - y h ). (1.2) 

i=l 
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This system determines deformations of the differential operator: 



V x :=^-Af(X), (1-3) 



n N- 



where 

n 

\F(\) ■= V + V 

. =1 A - a, 

that preserve its monodromy around the regular singular points {A = a.j}j = i 5 ... n and at A = oo. 
It was observed in [JMMS] (appendix 5) that, expressing the iVj's as 

Nt = GfFi, (1.5) 

where {F it Gi G M^ Xr } l=1) . ,. n are pairs of maximal rank kiXr rectangular matrices (ki < r), 
with {FiGj = Li G y/(fci)}i=i,... n constant matrices related to the monodromy of T>\ at 
{cti}i=i,...ni e Q- (1-1) ma y be expressed as a set of compatible nonautonomous Hamiltonian 
systems: 

dFi = {Fi,u)} (1.6a) 
dGi = {Gi,oj}. (1.66) 

Here d denotes the total differential with respect to the variables {ctj, y a } t=i,... n , the 1— form 

a — 1, . . . r 

to is defined as: 

n r 

w :=^2Hidai + ^2K a dy a , (1.7) 



i=l a=l 



with 



ffi := tr^) + y tT ( NiN j\ i=l,..., w (1.8a) 
:= V ai(Ni) aa + > =£ , a = l,...,r (1.86). 

^ Va-Vb 

bjta 

and the Poisson brackets in the space of (Fj,Gj)'s are defined to be such that the matrix 
elements of {Fj, Gi}i=i,... n are canonically conjugate: 

{(Fi) a . a , (Gj) M } = SijS a . b .Sab (1.9) 
i,j = l,...,n, a,b=l,...,r, a i? 6» = 1, . . . , fe». 
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The Frobenius integrability of the differential system (1.1) follows from the fact that the 
Hamiltonians {H^, K a } »=i,...„ Poisson commute. It was also noted in [JMMS] that, with 

a — 1 ,...r 

respect to the Poisson brackets (1.9), the matrices {Ni} i= x,...,n satisfy: 

{{Ni)abi {Nj)cd} = Sij {S bc {Ni) 

ad $ad 

(Ni) cb ), (1.10) 

which is the Lie Poisson bracket on the space (®f =1 gl(r))* dual to the direct sum of n copies 
of the Lie algebra gl(r) with itself. The 1-form u> is exact on the parameter space and may be 
interpreted as the logarithmic differential of the r-function, 

U = cflogT. (1.11) 

Numerous applications of such systems exist; in particular, in the calculation of correla- 
tion functions for integrable models in statistical mechanics and quantum field theory [JMMS, 
IIKS], in matrix models of two dimensional quantum gravity [M] and in the computation of 
level spacing distribution functions for random matrix ensembles [TW1, TW2], 

In the next section these systems will be examined within the context of loop algebras, 
using an approach originally developed for the autonomous case, involving isospectral flows, 
in [AHP, AHH1]. This is based on "dual" pairs of parameter dependent moment maps from 
symplectic vector spaces to two different loop algebras gl(r) and gl(N), where N = J27=i 
The nonautonomous Hamiltonian systems (1.6)-(1.8) will be generated by pulling back certain 
spectral invariants, viewed as polynomial functions on rational coadjoint orbits, under these 
moment maps, and identifying the parameters determining the maps with the deformation 
parameters of the system. This construction leads to a pair of "dual" first order matrix differ- 
ential operators with regular singular points at finite values of the spectral parameter, both 
of whose monodromy data are invariant under the deformations generated by these Hamil- 
tonian systems. In section 3, the generic systems so obtained will be reduced under various 
discrete and continuous Hamiltonian symmetry groups. A rank r = 2s generalization of the 
systems determining the correlation functions for an impenetrable Bose gas (or, equivalently, 
the generating function for the level spacing distribution functions for random matrix ensem- 
bles [TW1]) will be derived by reduction to the symplectic loop algebra sp(2s). The "dual" 
isomonodromy representations of the equations for the Painleve transcendents Py and Pyj 
will also be derived, and their Hamiltonian structure deduced through reductions under con- 
tinuous groups. A brief discussion of generalizations to systems with irregular singular points 
is given in section 4. 

2. Loop Algebra Moment Maps and Spectral Invariants 
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In [AHP, AHH1], an approach to the embedding of finite dimensional integrable sys- 
tems into rational coadjoint orbits of loop algebras was developed, based on a parametric 
family of equivariant moment maps J A '■ M — > gl( r ) + from the space M = {F, G G M Nxr } 
of pairs of N x r rectangular matrices, with canonical symplectic structure 

u = tr(dFAdG T ) (2.1) 

to the dual of the positive half of the loop algebra gl(r). The maps J a, which are parametrized 
by the choice of an N x N matrix A G M NxN with generalized eigenspaces of dimension 
{ki}i=i,...n, Y!i=i h = N and eigenvalues {a i } i = lj ... n , are defined by: 

J A :(F,G)\ — >G T (A-XI N )- 1 F (2.2) 

where In denotes the N x N identity matrix. The conventions here are such that all the 
eigenvalues {ai}i=i,... n are interior to a circle S 1 in the complex A-plane on which the loop 
algebra elements X(X) G gl(r) are defined. The two subalgebras gl(r) + ,gl(r)- consist of 
elements X + G gl(r)+, X- G gl{r)- that admit holomorphic extensions, respectively, to the 
interior and exterior regions, with X_(oo) = 0. The space gl(r) is identified as a dense 
subspace of its dual space gl(r)*, through the pairing 

<X U X 2 > = <f ti (X 1 (X)X 2 (X)) dX (2.3) 
Js 1 

X 1 egl(r)*, X 2 egl(r). 

This also gives identifications of the dual spaces gl(r)± with the opposite subalgebras gl(r) T . 
Taking the simplest case, when A is diagonal, the image of the moment map is 

AT. 

Afo(X) = G T (A - XI^-'F = V -— (2.4a) 

< * A — rv : 



i=l 



Ni := -Gf Fi, (2.46) 

where (Fi,Gi) are the k\ x r blocks in (F,G) corresponding to the eigenspaces of A with 
eigenvalues {(Xi}i=i,... n - The set of all Ao G gl(r)- having the pole structure given in eq. (2.4a) 
forms a Poisson subspace of o/(r)_, which we denote g^. The coadjoint action of the loop 
group Gl(r) + corresponding to the algebra gl(r) + , restricted to the subspace g^, is given by: 

Ad*(Gl(r) + ):gA^gA 

9i'=9{ai), i = l,...,n. 
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We see that could equally have been identified with the dual space (®™ =1 gl(r))* of the 
direct sum of n copies of gl(r) with itself, and the Lie Poisson bracket on gl(r)+ ~ gl(r)_: 

{/, g}W =<^W,dg]W >, (2.6) 

reduces on the Poisson subspace g^ to that for (®^ =1 y/(r))*, as given in eq. (1.10). 

In the approach developed in [AHP, AHH1], one studies commuting Hamiltonian flows 
on spaces of type gA (in general, rational Poisson subspaces involving higher order poles if 
the matrix A is nondiagonalizable) , generated by elements of the Poisson commuting spectral 
ring X\ of polynomials on gl(r)* invariant under the Ad*G/(r)-action (conjugation by loop 
group elements), restricted to the affine subspace Y + g^, where Y e gl(r) is some fixed 
r x r matrix. The pullback of such Hamiltonians under J a generates commuting flows in 
M that project to the quotient of M by the Hamiltonian action of the stability subgroup 
G A ■= Stab(A) C Gl(N). The Adler-Kostant-Symes (AKS) theorem then tells us that: 

(i) Any two elements of X\ Poisson commute (and hence, so do their pullbacks under the 
Poisson map J a)- 

(ii) Hamilton's equations for H e l\ have the Lax pair form: 



dt 
where 



f/A [(dH) + ,Af] = -[(dH)_,Af], (2.7) 



Af(\,t) := Y +Af (\,t), (2.8) 



with Ao £ gK r )- °f the form (2.4a), dH\j^ viewed as an element of (gl(r)*)* ~ gl(r), 
and the subscripts ± denoting projections to the subspace gl(r)±. 

The coefficients of the spectral curve of A/"(A), determined by the characteristic equation 

det(y + G T (A-XI N )- 1 -zl r ) =0, (2.9) 

are the generators of the ring X\ . 
In particular, choosing 

Y -diagd/!,...,^) (2.10) 

and defining {Hi e l\}i=i,... n by: 

Hi(M) :=^-.l tr((A/-(A)) 2 dA = tr(YA^) + V t AMA ( 2 .H) 
4ttz J x=a% cti - ctj 
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(where § x=a . denotes integration around a small loop containing only this pole), we see that 
these coincide with the H^s defined in eq. (1.8a). Thus, the Poisson commutativity of the 
Hi's follows from the AKS theorem. The Lax form of Hamilton's equations is: 



dt 

where 



dX :-[(dHi)-,W (2-12) 



(dHi). = -^-egl(r)-. (2.13) 
A - oti 



Evaluating residues at each A = a i: we see that this is equivalent to: 

dNj _ [N^N, 



8t{ OLj — OLi 



, j^h i,j = l,...n. (2.14a) 



= [Y + Y" '—m- (2-146) 

f * rv ■ — rv ■ 



Qti ^ OL^ 

If we now identify the flow parameters {ti}«=i, ■ n with the eigenvalues {on}i=i,... n , we 
obtain the nonautonomous Hamiltonian systems 

m l= [N t N^ f ^. hJ = 1 ^ ^ (215a) 

OCti OLj — OLi 

which are the OLi components of the system (1.1). Viewing the iVVs as functions on the 
fixed phase space M, eqs. (2.15a,b) are induced by the nonautonomous Hamiltonian sys- 
tems generated by the pullback of the H^s under the parameter dependent moment map J a- 
Eqs. (2.15a,b) are equivalent to replacing the Lax equations (2.12) by the system: 

which is just the condition of commutativity of the system of operators {V\, T>i\ i= i... n , where 
V\ is given by (1.3) and 

d d N- 

V t :=■£- + (dff,)-(A) = 7T- + t— (2.17) 

OCii OOLi A — OLi 

Remark. The system (2.16) could also be viewed as a Lax equation defined on the dual 
of the centrally extended loop algebra gl(r) A : in which the Ad* action is given by gauge 
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transformations rather than conjugation [RS]. The analogue of the spectral ring X\ is the 
ring of monodromy invariants, restricted to a suitable Poisson subspace with respect to a 
modified (.R-matrix) Lie Poisson bracket structure. This viewpoint will not be developed 
here, but is essential to deriving such systems through reductions of autonomous Hamiltonian 
systems of PDE's. 

The fact that the matrices {FiGf = Li e <7^(fti)}i=i,...n are constant under the defor- 
mations generated by the eqs. (2.15a,b), (2.16) follows from the fact that 

J Ga (F, G) := diag(F 1 Gf , . . . , F n G T n ) G gl(N) (2.18) 

is the moment map generating the Hamiltonian action of the stabilizer of A in Gl(N): 

N 

G a := JJ Gl(ki) = Stab(A) C Gl(N), (2.19) 

i=l 

this action being given by 

G A :M — > M 
K : (F, G) — > (KF, {K T )- X G) (2.20) 
K = diag(K 1 ,...,K n ) eG A , KiEGlih). 

The orbits of Ga are just the fibres of J a, so (J a, Jg a ) form a "dual pair" of moment maps 
[W]. Evidently, the pullback J\{H) is -invariant for all H e T\, and hence Jq a is constant 
under the Hi flows. 

So far, we have only considered the part of the system (1.1) relating to the parameters 
{<Xi}i=i,...n- What about the Hamiltonians {K a } a= i r that generate the y a components? As 
shown in [AHH1], besides J a there is, for each Y e gl{r), another moment map 

J Y :M^gl(N)X~gl(N)_ 
J Y : (F, G) i — ► -F(Y - zI r y x G T , (2.21) 

where z denotes the loop parameter for the loop algebra gl(N), whose elements are defined on 
a circle S 1 in the complex z-plane containing the eigenvalues of Y in its interior. The pairing 
identifying gl(N) as a dense subspace of gl(N)* is defined similarly to (2.3), for elements 
Xi E gl(N)*, X 2 € gl(N). The subalgebras gl(N)± are similarly defined with respect to this 
circle, and their dual spaces gl(N)± are identified analogously with gl(N) T . 

The moment map Jy is also "dual" to J a, but in a different sense than Jq a - one that 
is relevant for the construction of the remaining Hamiltonians {if a }a=i,...r- The space 
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may be identified with the quotient Poisson manifold M/Ga, with symplectic leaves given by 
the level sets of the symmetric invariants formed from each FiGj ', since these are the Casimir 
invariants on g^. Since the Hamiltonians in l\ are all also invariant under the action of the 
stabilizer Gy = Stab(Y) C Gl(r), where Gl(r) C Gl(r) is the subgoup of constant loops, 
we may also quotient by this action to obtain g^/Gy = M/(Gy x Ga)- Doing this in the 
opposite order, we may define gy C gl(N)+ ~ gl(N)- as the Poisson subspace consisting of 
elements of the form: 

M (z) = —F(Y - zI N )~ l G T = Y (2.22a) 

{M a ) ir .= F ia G ja , i,j = l,...,n, a = l,...r (2.226) 

(where, if the {y a }a=i,...r are distinct, the residue matrices are all of rank 1), and identify gy 
with M/Gy. Denoting by Xy the ring of Ad*-invariant polynomials on gl(N)*, restricted to 
the affine subspace —A + gy consisting of elements of the form 

M = -A + M , M egy, (2.23) 

the pullback of the ring Xy under the moment map Jy also gives a Poisson commuting ring 
whose elements are both Gy and Ga invariant, and hence project to Mj {Gy x Ga)- In fact, 
the two rings J A {X\) and Jy(Xy) coincide (cf. [AHH1]), because of the identity: 

det(A - AXv)det(y + G T (A - Mn^F - zl r ) = det(Y - zI r )det(A + F(Y - zI r y 1 G T - XI N ), 

(2.24) 

which shows that the spectral curve of Ai(z), defined by 

det(A + F(Y - zlr^G 7 - XI N ) = (2.25) 

and that of A/"(A), given by eq. (2.9), are birationally equivalent (after reducing out the trivial 
factors det(A — XI n) and det(Y — zl r )). 

Now, similarly to the definition of the elements {Hi e X^}j = i ... n , we may define {K a E 
1y}a=i,...r as: 

K a :=^-l tr(M(z)) 2 dz = -tv(AM a ) + V t AM^l (2 .26) 

*™Jz=v a Va-Vb 

To verify that these coincide with the K a 's defined in eq. (1.8b), we use eqs. (2.4a), (2.8), 
(2.22a) and (2.23) to express K a as: 

K a = — l dz — I dXXtr(M(z)(A- XIn)' 1 ) 2 (2.27a) 

47U J z=ya 2-Kl Jgi 

= -r—.i dXX— I dz [tr((y - zJ r )-W(A)) 2 - 2tr ((Y - 2/ r )"W(A))] .(2.27) 

4tT* J S 1 ^ % J z=y a 
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Evaluating residues at {y a } a =i,...r m z and at oo in A gives (1.8b). The Poisson commutativity 
of the K a } s again follows from the AKS theorem, and the commutativity with the H^s follows 
from the equality of the two rings J%^A = Jy^-Y- To compute the Lax form of the equations 
of motion generated by the K^s, we evaluate their differentials, viewing them as functions of 
jV(A) defined by eq. (2.27). Evaluating the z integral, this gives: 

dK M = XE a + ± E "^ E > + E ^ E " e 9 » (2.28) 

b ^ 1 Va Vb 

b^a 

where E a denotes the elementary diagonal r x r matrix with (aa) entry equal to 1 and zeroes 
elsewhere. Taking the projection to gl(r) + gives: 

(dK a ) + (X) = XE a + ± f f*W b + E b N iEa e gl(r) + , (2.29) 

and hence 

r 

^2(dK a ) + dy a = XY + Q, (2.30) 

where © is defined in eq. (1.2). By the AKS theorem, the autonomous form of the equations 
of motion is 

^L = [(dK a ) + ,Af], (2.31) 

while the nonautonomous version is 

jj£ = l(dK a ) + ,Af] + = [(dK a ) + ,Af] + E a . (2.32) 

Evaluating residues at A = ctj gives the equations 

dNj 
dr a 

for the autonomous case and 



[(dK a ) + ( ai ) t Ni], i = l,...,n (2.33) 



dy a 



[(dKj+ia&Ni], z = l,...,n (2.34) 



for the nonautonomous one. Eqs. (2.34) are just the y a components of the system (1.1). 
Eq. (2.32) is equivalent to the commutativity of the operators {T>\, X > o}o=i,...d where 

V:-.= -^--(dK a ) + (X), (2.35) 
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and implies that the monodromy of V\ is invariant under the y a deformations. In fact, it fol- 
lows from the AKS theorem that the complete system of operators {D\, Dj, X > a}i=i,... n ,o=i,...r 
commutes. 

Turning now to the dual system, it follows from the AKS theorem that the Lax form of 
the equations of motion induced by the K a : s on gl(N)-, viewed now as functions of M, in 
the autonomous case is 

dM 

^ = -[(dK a )_,M], (2.36) 

where 

(dK a )_(z) = G gl(N)_. (2.37) 

z-y a 

(Note that the differential dK a entering in eqs. (2.36), (2.37) and below has a different signif- 
icance from that appearing in eqs. (2.28)-(2.35).) Evaluating residues at z = y a shows that 
(2.36) is equivalent to the system 

dM b [M b ,M a ] 

— — = -, b f: a, a,o=l,...r, (2.38a) 

OTa Vb-Va 

<~> I a b=1 IJa Ub 

b^a 

Identifying the flow parameters {r a } a= i j ... r now with the eigenvalues {y a }a=i,...r °f Y gives 
the nonautonomous Hamiltonian system 

dM b [M b , M a ] 

— — = -, b f: a, a,b=l,...r, (2.39a) 

oy a Vb - y a 

Equations (2.39a,b), (2.40) are equivalent to the commutativity of the system of operators 
{T>z,1>a}a=i,...r defined by: 

V z := — -M(z), (2.41a) 
B . := £ + (,«.)_(,) = £ + _^_, a= (2Alb] 

Thus the monodromy of the "dual" operator V z: is also preserved under the y a deformations. 



or, equivalent ly, 
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Finally, by similar computations to the above, with (A, A, AT) and (Y, z, M) interchanged, 
we can express the H^s as functions of M: 

Hi=^-l dzz^- <f dX [tr((A - XI n)' 1 M{z)f - 2tr ((A - XIn^M^))] . (2.42) 
4iti Jgi 2tti J \ =ct . 

Evaluating residues at {aj}j=i,... n in X and at oo in z gives 

H t = - V(M )« + V E=iWi J -)(ELiW)ii) > 1 = 1 ,..., n . (2 .43) 

To compute the Lax equations for .M(z) generated by the i^i's, we evaluate their differentials 
when viewed as functions on gl(N)- defined by eq. (2.42). Evaluating the A integral, this 
gives: 

mW = -zE, + ± E > M{Z)E > + E > M{Z)E > 6 fox), (2.44) 

where E'i now denotes the elementary diagonal N x N matrix with (ii) entry equal to 1 and 
zeroes elsewhere. Taking the projection to gl(N) + gives: 

(dH iM z) = _ zEi + j-± E iM^l±ElM^i 6 S 1 (JV)+ , (2.45) 

J=l a=l 1 J 

and hence 

n 

^2(dH l ) + (z)da l = -zdA + (2.46) 
i=i 

where 

r 

$ y = (1 - <y (J](M a )^log(a, - a,). (2.47) 

a=l 

By the AKS theorem, the autonomous form of the equations of motion is 

dM 



while the nonautonomous system is 



[(dHi)+,M] (2.48) 



9M = [(dHi)+,M] + d{ t Hi)+ = [(dHi)+,M] - E t . (2.49) 



dai dz 
Evaluating residues at A = a.i gives the equations 



^ = [(dH i ) + (y a ),M a ] (2.50) 
11 



for the autonomous case and 

^ = [(dH i ) + (y a ),M a ] (2.51) 

for the nonautonomous one. Eq. (2.49) is equivalent to the commutativity of the operators 
{V z ,V*} i=1 _ n , where 

V*:=-^--(dH i ) + (z), (2.52) 

and implies that the monodromy of V z is invariant under the ai deformations. Again, it fol- 
lows from the AKS theorem that the complete system of operators {V z , V a , £ > *} a =i,...r,i=i,...n 
commutes. 

Thus, at the level of the reduced spaces gA/Gy ~ gy/GA, we have two equivalent 
"dual" isomonodromy representations of the Hamiltonian systems generated by the spectral 
invariants {Hi, -Ka}i=i,...,n,a=i,...,r - systems (2.15a,b), (2.24) intheA/"(A) G gA representation 
and (2.39a,b), (2.51) in the M.{z) G gy representation. 

3. Reductions 

The general scheme of [AHP, AHH1] may combined with continuous or discrete Hamil- 
tonian symmetry reductions to deduce systems corresponding to subalgebras of gl(r) and gl(N) 
or, more generally, to invariant submanifolds. In particular, the Marsden-Weinstein reduction 
method may be applied to the symmetry groups Ga and Gy, or to other invariants of the 
system. 

The discrete reduction method (see [AHP], [HHM] for further details) may be sum- 
marized as follows. Suppose a r * : gl( r ) + ~ * 9K r )+ ^ s a Lie algebra homomorphism that is 
semisimple, of finite order and induced by the group homomorphism a r : GZ(r)+ — > Gl(r) + . 
Let a* : gl(r)- — * gl(r)- denote the dual map, which is a Poisson homomorphism, and let 
h a C gl(r) + , h* C gl{r)_ denote the subspaces consisting of the fixed point sets under a r * 
and cr*, respectively. (These are naturally dual to each other, since h* may be identified with 
the annihilator of the complement of h a under the decomposition of gl(r)+ into eigenspaces 
of cr r *)- Then h a C gl(r)+ is a subalgebra, and its dual space h* a has the corresponding Lie 
Poisson structure. Suppose there also exists a finite order symplectomorphism o~m : M — > M 
such that the moment map J a satisfies the intertwining property: 

J A o (TM = &* o J A- (3.1) 

The fixed point set M a C M is, generally, a symplectic submanifold, invariant under the flows 
generated by o"M-invariant Hamiltonians. The restriction Ja\m„ '■= J Aa takes its values in 
h* a , defining a Poisson map: 

J A , : M a — > hi, (3.2) 
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which is the equivarient moment map generating the action of the subgroup H a C Gl(r)+ 
consisting of the fixed points under 5> : Gl(r)+ — > Gl(r) + . Such reductions, when applied 
to the spectral invariants on Y + gA generate systems satisfying the criteria of the AKS 
theorem, provided the matrices Y and A are appropriately chosen to be compatible with the 
reduction. The same procedure may be applied to the dual systems on —A + gy if a similar 
homomorphism cx/v* : gl(N) + — » gl(N) + exists, satisfying the intertwining property: 

Jy ° o~m = &n ° Jy- (3-3) 



The corresponding moment map, 

J Ya : M a — k*, (3.4) 

obtained by restriction Jy a := Jy\m„ takes its values in the Poisson subspace k* C gl(N)_ 
consisting of the fixed point set under the dual map a* N : gl(N)- — > gl(N)-, and k a C gl(N)+ 
is the corresponding subalgebra consisting of fixed points under ctn*- 

The following examples illustrate both the discrete and continuous reduction procedures. 

3a. Symplectic Reduction (discrete) 

Let r = 2s and define a*,a* N and o~m by: 

a* : X(X) — ► JX T (X)J, (3.5a) 
a* N :Z(z)^e(-z), (3.56) 
om : (F, G) — ► J(GJ, -FJ), (3.5c) 

Xegl(r)_, Zegl(N)_, (F, G) e M 



J = I T ° ) e M 2sx2s . (3.6) 



where 

T=( 

, -Is 

Then M a C M consists of pairs (F, G) of the form: 

F=^=(Q,P), G=^=(P,-Q) (3.7) 

where Q,P £ M Nxs are N x s matrices. The blocks {(Fi, Gi)}i=i,... n corresponding to the 
eigenvalues {0^)^=1,. .. n are similarly of the form 

F i = -j=(Q i ,P i ), G i =-j=(P i ,-Q i ), (3.8) 
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where Qi,Pi G M kiXs . The reduced symplectic form on M a is 

u = tr (dQ A dP T ). (3.9) 

The subalgebra h a C gl(r)+ is just the positive part sp(2s) + of the symplectic loop algebra 
sp(2s), and the dual space h* is similarly identified with sp(2s)_. The subalgbra k a C gl(N) + 

~(2) ~(2) 

is the positive part (7/ (-/V)+ of the "twisted" loop algebra gl (N), with dual space k* ~ 

~(2) ~ 

(7/ (-/V)_. The image of the moment map has the form 

-PfQ* -pTp. 

AT (A) = J^(g,P) = £ V ^ Qf A 7 ■ (3.10) 

i=i 1 

In order that the pullback of the elements of the ring under J a be o"M-invariant, and that 
the intertwining property (3.3) for the moment map Jy be satisfied, the matrix Y must be in 
sp(2s). For diagonal Y, this means 

y = il °A (3-11) 



-y 

V = diag(yi, ... y a ) G gl(s). 
The image of Jyv then has the block form 

A,„ W = *, W ,^£(^-£|), (3,2) 

where {q a } a =i,...r and {p a }a=i,...r denote the ath columns of Q and P, respectively. 
The Hamiltonians {#i}i=i,... n reduce in this case to: 

^-M^,),^" 1 ^-^^-^' 1 (3.13) 
and generate the equations of motion: 



j ^ i (3.14a) 



dQ, _ PjQjQj - QjP[Qj 
dai 2{a.i — ctj) 

^ = AO^-^ , (3 , 46) 
oati 2(ai — aj) 

9Q, n , ^ QiPfQj - PQJQ, .... . 

^- = -Q^ + E 2(a< _ aj) (3-14.) 

-K—=P l y + 2^ ^77 r — • (3.14d) 
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The Hamiltonians {K a = — if a +s}a=i,...s reduce to: 



K. = \?lA qa + \ ± - 1 ± *lwfa£) (3 . 1B) 



and generate the equations 



b^a 



^4^ + 1^, b^a (3.16a) 
oy b Ay b -y a 4y b + y a 

9|/6 4 y b - y a 4 y 6 + y a 

^ = - i V + 1 y (3 .i6d) 

b^a b^a 



The particular case s = l,{ki = l}i=i,... n of (3.13), (3.14a-d) reduces, up to a simple 
change of basis, to the system of Theorem 7.5, [JMMS]. The corresponding r-function gives 
the n-particle correlation function for an impenetrable Bose gas or, equivalently, the level 
spacing distribution function for a set of random matrices having no eigenvalues in the intervals 
{[cK2i-i, Oi2i]}i=i,...m, n = 2m, in the scaling limit [TW1] . 

Hamiltonian Structure of Painleve Equations 

The following two examples show how the Painleve transcendents Py and Pyi may be 
derived from the generic systems (2.15a,b), (2.34) or (2.39a,b), (2.51) through Hamiltonian re- 
duction under continuous symmetry groups. Our derivation will be guided by the formulation 
of Painleve transcendents as monodromy preserving deformation equations given in [JM], but 
the emphasis here will be on the loop algebra content, the Hamiltonian reductions and the 
associated "dual" systems. For previous work on the Hamiltonian structure of the Painleve 
transcendents, see [Ok] and references therein. 

3b. Painleve V 

Choose N = 2, r = 2, and 
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Then F and G are 2x2 matrices 

*-(«)• G =(S)' < 3 - 18 > 

with row vectors {Fj = (Fn F i2 ), Gi = (Gn Gi 2 )}i = i, 2 . The stabilizer Ga C G7(2) of A is 
the diagonal subgroup generated by the moment map 

J Ga (F,G) = {F 1 Gl,F 2 Gl) := ( m ,^ 2 ). (3.19) 

Fixing a level set, we parametrize the quotient under this abelian Hamiltonian group action 
by choosing the symplectic section Ma C M defined (on a suitable connected component) by 

i f X! yi-%\ n_ i fyi + ^ -xi 

M2 



V2\xi »-gy G "v^U + S (3 ' 20) 

The reduced manifold M re( j = Jg* (/zi, h 2 )/Ga is identified with R 2 x R 2 minus the coordinate 

axes {xi = 0, x 2 = 0}, quotiented by the group of reflections in these axes. The reduced 
symplectic form is 

2 

^red = dxi A e%. (3.21) 

The image of the reduced moment map J a '■ Af re d — > <7^(2)_ translated by Y - is 
JV(A) = y + J A (F,G) 



* 



M 2 \ / M 2 



"2 

X 2 XlJ/1— Ml / _|_ V »2 



(3.22) 



-ty 2A 2(A-1) 

The stabilizer GV C G7(2) of Y is the diagonal subgroup, acting by conjugation on Ai(X), 
which corresponds to scaling transformations 

(xi,x 2 , yi, V2) — > (e r xi, e T x 2 , e" T yi, e~ r y 2 ), (3.23) 

and is generated by 

a-=-(xiy 1 + x 2 y2)- (3.24) 
(The trace part of gl(2) acts trivially, since it coincides with the Casimir + /j, 2 .) 
The Gy-invariant Hamiltonians Hi , H 2 e may be expressed 

Hi =-*- / tr(W(A)) 2 dA = tr(YATi) - tr(iViiV 2 ) = -tH - - a 2 - 2at (3.25a) 

— / tr(A/"(A)) 2 dA = tr(Y/V 2 ) + tv(NiN 2 ) = «f + + a 2 . (3.256) 

7™ 7A=1 2 



H 2 . 

4"7U 
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where 



H = -^ + x2 ^ + y 2 ^ + li ~ X2V2 ' (3 ' 26) 



The dual system is determined by the moment map Jy '■ M — > (7/(2) defined by eq. (2.21) 
which, when restricted to the symplectic section Ma defined by eq. (3.20), gives 



M(z) = - A - F(Y - zI 2 )~ L G 



xiyi+m x 1 y 2 +^2^ \ / -xiyi+m -a^J/i+Miff- 



\) ' 2Tr^i + w+t) ' ( 

Here, the quantities a and \i\ + fj, 2 are interpreted as Casimir invariants, whereas fix and fi 2 
individually are conserved quantities because they belong to the spectral ring Zy. In terms 
of the dual system, we may express the Hamiltonians K l7 K 2 G ly as: 

K x = ^— I tr(M(z)) 2 dz = -tr(AMi) + — tr(MiM 2 ) 
4ttz J x=t 2t 

= E + d±A_^ (3.28a) 
2 8t 2 v y 



^2 = ~t~~t / tr(_M(z)) 2 d,2 = -tr{AM 2 ) - ^-tr(M 1 M 2 ) 
4ttz J x=t 2t 

H nl + n 2 2 [i 2 



2 8* 



(3.286) 



The relations 



H x = -H 2 - 2at (3.29a) 



(fix - (to) 



2 



= -2tKx + — 4 ^ ; - a 2 - 2at - fy/ 2 (3.296) 
= 2tK 2 + - a 2 - 2at + t/i 2 (3.29c) 

can also be derived from the identity (following from eq. (2.24)): 

z 2 - ztrAf(\) + \ ((trAT(A)) 2 - trA/" 2 (A)) _ A 2 + XtrM(z) + | ((trA^)) 2 - tr_M 2 (^)) 

z^t 2 ~ A(A-l) ' 

(3.30) 

Integrating both sides around contours in the z-plane containing either the pole at z = t or 
the one at z = — t and contours in the A-plane containing either the pole at A = or the one 
at 1, we obtain eqs. (3.29a-c). 
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Viewing Ki and K 2 as functions of M we have, from eq. (2.29) 

<«>>-(s SHU* ^ 
^-^s 



and hence the monodromy preserving deformation equation for A/"(A) generated by K\ — K 2 

d_ 



is given by the commutativity of the operators V\ = —J\f(\) and T>\, where M{X) is given 
by eq. (3.22), and 



n = ^~( X n ->'-^^\. (3.32) 



dt \0 -X J 2t \xi+x 2 2 
Viewing K\ and K 2 instead as functions of M, from (2.37) we have 

(d Kl )_ = ^— (3.33a) 

— aiyi+Mi -x 2 yi+fj-i^- 
(dK 2 )_ = ^_ . (3.336) 

To obtain the corresponding dual monodromy preserving deformation equation for Ai(z), we 
cannot simply restrict eqs. (2.39a,b), (2.40) to the submanifold Ma C M. The image of Ma 
under Jy is 3-dimensional, since points related by the scaling transformation (3.23) have the 
same image. Unlike Ja\m a , the restriction Jy\m a is not a Poisson map, since Ga does not leave 
Jy invariant, but acts by conjugation on the image Jy(F, G). Therefore Jy does not project 
to define a map on the quotient space M/Ga- However, the reduced systems generated by 
the Ga and GV-invariant Hamiltonians H±, H 2 , K\ and K 2 are determined by the projection 
of their Lax equations to the quotient manifold M/Ga- The projected (nonautonomous) 
Hamiltonian vector field determined by K\ — K 2 has a unique lift that preserves the section 
Ma, obtained by adding a "vertical" term 

dK v := 1 2 3.34 

to the factor [dK\)- — (dK 2 )- entering in the Lax equation for M{z). Eq. (3.34) is obtained 
by noting that, apart from the conditions that the diagonal terms [i\,[i 2 in M\ + M 2 be 
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conserved, and the Casimirs detMi , detM 2 vanish, all of which are automatically satisfied by 
the Lax system, the only remaining condition defining the image Jy(Ma) is: 

det (M(z) - M T (-z)) = 0, (3.35a) 

or, equivalent ly, 

det (Mi + Ml) = 0. (3.356) 

Up to multiples of the identity matrix, the unique element in the diagonal subalgebra (cor- 
responding to Ga) which, when added to (dKi)- — (dK 2 )-, preserves this condition is dK v . 
Thus, the correct dual deformation operator T> t , whose commutativity with V z — J| — Ai(z), 
gives the Hamiltonian system generated by K\ — K 2: is: 



\ +Mlf^ J \-Xiy2+H2^ -X2V2+V2 



dt 2{z-t) 2(z + t) 

+ 



1 / Mia: 2 \ r, 

M-zr—^) 
x i 

2 



(3.36) 



1 / V2 X \ \ 
2 



To reduce the system under the Gy action generated by a, we first introduce the "spec- 
tral Darboux coordinates" (-U, vS) (cf. [AHH2]) defined by 

where 

w = x\ + x\ (3.38a) 

u = -^L (3.386) 

In terms of these, the symplectic form (3.21) becomes 

<^red = dlog-u; A da + du A dv, (3.39) 

where 

U^ + ^uA (340) 



2 \ u u — 1 r 

is the momentum conjugate to u and a is the invariant defined in eq. (3.24). The operators 
V\,V Z ,T>1 and V t may be expressed in terms of these coordinates by substituting 

x\ = uw, x 2 = (1 — u)w, Xiyi = 2u(a — uv + v), £22/2 = 2(1 — u)(a — uv ) 

9 a H / \9 9 a (1 ^0 / \9 (^'^"0 

yi = 4— (a - uv + v) z , y 2 = 4- -(a-uvY 

w w 
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in eqs. (3.22), (3.27), (3.32) and (3.36). 

Choosing a level set a = a , the symplectic form (3.21) reduces to 

Wred|a=a = du A dv i (3A2) 

so (u, v) provide canonical coordinates on the the reduced space obtained by quotienting by 
the Gy-flow. From eqs. (3.22), (3.25a,b) it follows that we may write 

From (3.37) and (3.40) it follows that 

l t ^ (u) = {v - t f + l(^ + J^)\ ,3.44) 

Evaluating the integral 

A(A - l)trA/" 2 (A) 1 



1 / A(A - l)trA/ 2 (A) 1 . 1U ,. 2 . , 

* — - 1 —d\ = -u(u - l)trM 2 (u) 

47tiJ x=u X-u 2 



* 2 u( u - 1) + u(# i + #2) - Hi + t ~ 1 \ " + 



2u 2(u-l) 

(3.45) 

and using eqs. (3.29a-c), (3.44), (3.45) gives 

Kl - K> = - 2„ t ) + 2au + i| - ^ 4 - 2, (3.46) 

(Note that the simple canonical change of coordinates v — > u + ^ + 2 (^_i) transforms this to 
polynomial form; cf. [Ok]). The reduced equations of motion generated by the Hamiltonian 
Ki — K 2 are therefore 

du = 2u{u-l) {v _ l) {SA7a) 

Qjty t 

— = (v 2 -2^) + -^ ; P2 -2a. 3.476 

dt t v y 4w 2 t 4(w-l) 2 t v ; 

Eliminating v by taking second derivatives gives: 



d 2 u ( 1 1 \ / du\ 2 1 du u u — 1 



dt 2 \2u 2{u-l))\dt) t dt t 2 {u-l) t 2 u 
_ 1 < u ~ 1 ) -Su{u-l){2u-l) 



(3.48) 
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where 

2 2 

a=^, /3 = -y, 7 = 4a + 2, 5 = 2. (3.49) 

which is one of the equivalent forms of Py. The more usual form is obtained by transforming 
to the new variable 

w = — ^— . (3.50) 
it — 1 

5c. Painleve VI 

We take N = 3, r = 2, A = diag(0, 1, t), Y = 0, so F, G are 3 x 2 matrices 

G=lG 2 l, (3.51) 

with rows {F; = (Fji F i2 ), Gj = (Gn Gi 2 )}i=i,...3. If t ^ 0, 1 the eigenvalues of A are distinct 
and the stabilizer Ga C GZ(3) is the diagonal subgroup generated by the moment map 

Jg a = (//i,// 2 ,// 3 ) = (^G^FaG^FaGi 1 ). (3.52) 

Fixing a level set, we parametrize the quotient under this abelian Hamiltonian group action 
by choosing the symplectic section Ma C M defined (on a suitable connected component) by 

(3.53) 

/ \ y + M / 

x 3 / \ y,i x 3 6 / 

(The choice of signs is made such that subsequent reductions be at nonsingular points in the 
real case). The reduced symplectic form is 

3 

^red = ^2dxi A dyi, (3.54) 
i=i 

and the reduced manifold M red = Jq a /x 2 , 1^3) /G a is identified with R 3 x R 3 minus the 

coordinate planes {x\ = 0, x-i = 0, £3 = 0}, quotiented by the group of reflections in these 

planes. The image of the reduced moment map J a '■ M re d — > (2) _ is 

2 2 
2 , **i \ / 2 ^9 

x\ X\V\—[1.\ J \ x\ X2V2-V2 





Af(\) = Afo(X) = J A (F, G) =^ 1 - iyi " w + 



2A 2(A-1) 



2 ^1 

-^3 2/3-M3 y 3 — f 

X 3 
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Choose the Hamiltonian 



H = H 3 I tr(A/"(A)) 2 dA 

4 ^ Jx=t 



tr(JViJV 3 ) tr(JViJV 3 ) 



t- 1 



1 

it 



2 3v O 



+ 



Let 



4(f-l) 
1 3 



{X 2 y3+X 3 y2) - A*|— | - //§-§ + 2/12//3 



6 =^ + ^-I/3 2 )-^ 
_ 1 / 2 , 2 2\ 



2 2 
A*2 , A*3 



2Xg 



(3.56) 

(3.57a) 

(3.576) 
(3.57c) 



(3.58) 
(3.59) 



These are the generators of the constant SI (2) conjugation action 

g:M (\)-^gAf (\)g-\ 

and satisfy 

{a, 6} = b, {c, a} = c, {6, c} = —2a. 

(In this case, since Y = 0, GV = G7(2), but the trace term acts trivially.) The Hamiltonian 
(3.56) is invariant under this S7(2)-action 

{a, H} = {b, H} = {c, H} = 0, (3.60) 

since the elements of the spectral ring X\ are Gy invariant. 

The monodromy preserving deformations generated by H are then determined by the 

JL 



commutativity of the operators T>\ = -j^ — M(\) and V t , with A/"(A) given by eq. (3.55) and 



d_ 
dt 



+ 



-X3V3-/J-3 



2 ^3 
3 

^32/3-^3 



(3.61) 



They also preserve the monodromy of the "dual" operator V z = -§■ — M.(z), where M(z) is 



2(X-t) 
il" opera 

determined by restricting the moment map Jy '■ M re a — > fl^(3)_ to the submanifold C M 
M (z)=-A + M (z) = -A + J Y (F, G) = 




1 



1 



2mi 



t 



xiy2-x 2 yi + ^f- + ^- 



2M2 



M1!C3 _ ^3^1 



^3 2/2+^2 2/3 + 



^2^3 _ ^3^2 



1 Ml ^3 1 
KlJ/3 +^32/1 j^H — ^ 

x 2 2/3+-32/2-^+^if 1(3.62) 

2M3 
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In this gl(3)~ representation, the quantities 

tr(FG T ) = m + ii2 + ^3 
tr(FG T ) 2 = tr(G T F) 2 = 2(a 2 - be) + ±(£ 
det(FG T ) = 



i=i 



(3.63a) 
(3.636) 
(3.63c) 



are the Casimir invariants, while the individual elements ^1,^2,^3 are not Casimirs, but 
generators of the stabilizer Ga C Gl(S) of A, and hence constants of motion. Thus, what 
appeared before as Casimirs on gl(2)_ become elements of the spectral ring Zy, while the 
element of Z\ given by eq. (3.63b) becomes a Casimir on gl(3)-. 

Viewing H3 now as a function of A4, we have, from eq. (2.45 ) 



(dH 3 ). 




1 

+ 2 





j(x 3 yi+x iy3 + ^3--^^) T r T (x 3 y 2 +x 2 y 3 + 



^ T (x 2 y 3 +x 3 y 2 -^ + ^) 



x 2 x 3 I 



(3.64) 

The dual monodromy preserving representation is therefore given by the commutativity of the 
operators V z — J| — M(z), and T>\, with M(z) given by eq. (3.62), and 

V* = ^-(dH 3 )+ + dH v 




( 



1 

2 



1/ M l x 3\ 
f(M3 -3T 2 -) 



1 1 ti i* 3 \ 



jh^y2+x 2y3 + ^-^ ) I( Ml -^i) +I A_( M5 



Here 



dH v = 



V 







"1 




2(t=T)^ ^) 





1 \ , 1 \ 

A(^!J i -Mi)+2(Fri)(^-M2)/ 



(3.65) 



(3.66) 



is the element of the diagonal subalgebra (corresponding to Ga) that must be added in order 
that the lift of the G^-reduced Hamiltonian vector field on M re d be tangential to Ma- 
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To obtain the <S7(2)-reduced system, first choose the level set 

b = c = 0, (3.67) 
and again define the "spectral Darboux coordinates" (u, w) by 

^M^+A-A)-^ (3 - 68) 

on this level set, where 



w = (1 +t)x\ + tx 2 2 - x\ (3.69a) 

u=-^- (3.6%) 
w 

a (A) := A(A - 1)(A - t). (3.69c) 
Thus (u, w) are elliptic-hyperbolic coordinates on the cone 

x\ + x\ = x\. (3.70) 
On the invariant manifold defined by b = c = 0, the symplectic form reduces to 

^red I (6,c)=(o,o) = d(logw) A da + du A dv, (3.71) 

where 

= 1 / gy/i - (li X 2 y 2 ~ ^2 £3?/3 - /o 

U 2 V u u - 1 u-t ) { ' ) 

is the momentum coordinate conjugate to u. (Note that the slight difference between this 
choice and that of eq. (3.40) results in a polynomial form for the Hamiltonian.) 

Restricting to the level set a = ao, we have 

^red|(a,6,c) = (a ,0,0) = du A dv, (3.73) 

which is the reduction of the symplectic form under the S7(2)-action generated by (a,b,c). 
The coordinates (u, v) project to the quotient under the action of the stability group of the 
image (ao,0, 0) of this sl(2) moment map, since they satisfy {u, a} = {v,a} = on the level 
set (a, 6, c) = (ao, 0, 0). 



To compute the Hamiltonian in terms of the reduced coordinates, we write 

l tr(m)f) = ™ + J± 2 + _d_ + (3 . 74) 
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where 



and evaluate 



\i=i j i=i 



-LI «(WA)) 2 dA= i ^ 2 

4ttz / A=tt A - u 2 

p.p. l)(u- t) {4u(u-t) fi 2 3 u(u-l) 

= Po + t'lU H h -— + —; — , 

where the integral is around a circle containing only the pole at A = u. Since 

^ + ^_ = 0, (3.77) 

u it — 1 u — t 



we have 

-■ 

2" v/ ' I « ' n-1 ' n-ti ' 2 I « ' u-1 ' u-t 



and hence, 



Pq =u(u — l)(u — t)v 2 + V (fJ>l(u — l)(u — t) + fJ*2u(u — t) + fJ<3u(u — 1)) 
+ Hl^2{u ~t)+ ^2^U + /Ui / U 3 (w - 1) - PlU. 



(3.79) 



From (3.56), (3.79) 

Pp + Pit 
t(t-l) 



H = 



[u(u — l)(li — t)v 2 + V (Hl(u — l)(li — t) + H2U{U — t) + ^3U(U — 1)) 



*(*-!)' 

+ ^i^iiu -t) + IX2IX3U + musiu -l) + (t- u)Pi]. (3.80) 
The reduced form of Hamilton's equations are therefore 

[2u(u — l)(it — t)v + (j,i(u — l)(it — t) + {i2u{u — t) + fJ,su(u — 1)] (3.81a) 



dt t(t - 1) 

+ v (/zi(2« - £ - 1) + jU 2 (2it -t) + /U 3 (2w - 1)) + A*iA*2 + + A*iA*3 - Pi]- 

(3.816) 
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Upon elimination of v, this gives Pyi- 



d 2 u 1 (\ I 1 \ fdu\ (\ 1 1 \ du 



dt 2 2 \u u — 1 u — t J \ dt J \t t — 1 u — t J dt 



u(u-l)(u-t) ( n t t-1 — 1) i 

+ %-V ( a+ ^ +7 (^TF + ^(biF)- (3 ' 82) 



where 

« = 2a 2 , 0=-^/?, 7=^/4 <*=-^3- (3-83) 



4. Generalizations 

The approach developed in [AHP, AHH1] is equally valid for matrices A, Y that are 
nondiagonalizable, giving rise to isospectral deformations of matrices of the form 

i=i i i= i v v 

r r " M 
o=lm„ = l^ yaJ 

A straightforward generalization of the moment map construction may also be made, yielding 
the more general forms 

l =0 i=l h = l {A ai > 

r r r a 



m =0 a=lm a =\ V yQj 



Such A/"(A), .M(z) may be viewed as elements, respectively, of subspaces C gl(r)*, g Y c 
gl(N)* defined by the rational structure appearing in eqs. (4.2a,b). These are Poisson sub- 
spaces with respect to the Lie Poisson bracket on gl(r)* (resp. gl(N)*) corresponding to the 
Lie bracket: 

[X,Y] R := ±[RX,Y] + ^[X,RY], (4.3) 

where 

R:=P+-P- (4.4) 

is the classical R— matrix given by the difference of the projection operators P± to the sub- 
algebras gl(r)± (resp. gl(N)±). The .R-matrix version of the AKS theorem [S] again implies 
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that results of the type (i) and (ii) (eq. (2.7)) hold for the autonomous systems generated by 
elements of the ring X\ (resp Xy) obtained by restriction of the ring of Ad*-invariant poly- 
nomials on gl(r)* (resp. gl(N)* to (resp. gy). Isomonodromic deformations of operators 
of the type T>\ = -J^ — jV(A), where jV(A) is of the form (4.2a), were the subject of the series 
of papers [JMU], [JM]. They are required, in particular, for the isomonodromy formulation 
of the remaining Painleve transcendent equations Pi - Piy ([JM]) and for the Hamiltonian 
dynamics governing the level-spacing distribution functions in random matrix models at the 
"edge of the spectrum" [TW2]. A brief discussion of the latter from the loop algebra view- 
point is given in [HTW] . The Hamiltonian formulation of more general systems of monodromy 
preserving deformations of operators with irregular singular points of arbitrary order within 
the framework of spectral invariants on loop algebras will be addressed in a subsequent work. 
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